PROCEEDINGS 

OF SCIENCE 



(N 

O 

(N 

C 
3 



(N 

S 1 



O 
=5 



> 

oo 

O 

(N 



Three-nucleon contact interaction at the 
next-to-leading order 



L. Girlanda 

Dipartimento di Matematica e Fisica "E. De Giorgi" 
Universita del Salento, and 
INFN Sezione di Lecce, 1-73100 Lecce, Italy 
E-mail: [girlanda@le . infn . it 



A. Kievsky 

INFN Sezione di Pisa, 1-56127 Pisa, Italy 



E-mail: kievsky@pi . infn . it 



M. Viviani 

INFN Sezione di Pisa, 1-56127 Pisa, Italy 



E-mail: viviani@pi.infn.it 



A satisfactory description of bound and scattering states of the three-nucleon (3N) system is still 
lacking, contrary to what happens in the two-nucleon case. In the framework of chiral effective 
theory, it is possible that a realistic 3N interaction will require the inclusion of subleading con- 
tact terms, which are unconstrained by chiral symmetry. We construct the minimal 3N contact 
Lagrangian imposing all constraints from the discrete symmetries, Fierz identities and Poincare 
covariance, and show that it consists of 10 independent operators. By integrating out the pions 
from the effective theory we derive the pion-induced contributions to the corresponding low- 
energy constants (LECs), and find no clear signal of convergence of the pionful theory between 
N2LO and N3LO. 
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1. Introduction 

One of the central tasks in nuclear physics is the determination of the interaction among nu- 
cleons. Traditionally, the two-nucleon (2A0 forces have been modeled phenomenologically. By 
now, these models reach a remarkable degree of accuracy, providing an excellent fit to the very 
large body of experimental data. More recently, so called chiral potentials have been developed, 
which establish a solid link between the nuclear interaction and the underlying strong interaction 
theory [|L ||, B, 0, B). They are derived in the framework of an effective field theory, incorporating 
all symmetry properties of QCD, and organized as a systematic expansion in powers of momenta 
and quark masses. This progress, together with the tremendous improvement in ab-initio numerical 
calculations of bound and scattering states of few-nucleon systems [g], make the latter a privileged 
theoretical laboratory to test our understanding of the nuclear interaction. 

Presently the 2N potential has been worked out up to the next-to-next-to-next-to leading order 
(N3LO) in the low-energy expansion [[| 0], yielding a % 2 per datum very close to one. In the 
three-nucleon (3A0 sector, however, the situation is not as satisfactory. A 3N interaction starts 
to show up at the N2LO [^, |Kj in the effective theoiy. From the phenomenological point of 
view, its contribution is necessary in order to reproduce the 3- and 4-nucleon binding energies. 
Nevertheless, all currently available 3N potentials fail to simultaneously describe also the neutron- 
deuteron doublet scattering length jTTj]. As a matter of fact, several low-energy data in the nucleon- 
deuteron system remain still unexplained, the most prominent one being the so-called A y puzzle. 

It is interesting to notice that the discrepancies between theory and experiment in the 3iV 
system arise especially at veiy low energy, and disappear as the energy increases. At such low 
energies, the nuclear interactions are effectively point-like. We therefore propose to refine that part 
of the interaction which stems from purely nucleonic vertices. Indeed, at very low energies, even 
the pions can be integrated out of the theory, giving rise to the "pionless" effective theoiy. The 
inclusion of contact interactions, which are unconstrained by chiral symmetry, could also provide 
the necessary flexibility to describe the 3N interaction. In this respect it is fitting to recall that all 
adopted 3N interaction models only contain a few free parameters, contrary to what happens for 
the 2N interaction, which is parametrized by more than 20 adjustable parameters. For example, 
the chiral 3N force at N2LO only contains 2 adjustable LECs, and its extension to N3LO does not 



involve any new LEC [J12J, [130. 

Invariance under parity dictates that subleading 3N contact terms contain two (covariant) 
derivatives of nucleon fields. In the pionful version of the effective theory such terms would there- 
fore contribute at N4LO, and are therefore beyond the accuracy of present day models, whereas 
in the pionless theory they represent the first corrections to the 3N interaction. It may well be 
that, at very low energy, the latter expansion scheme is more effective. Our task is twofold: first, 
classify all possible terms of this kind, second, fit the accompanying LECs to 3N data, in asso- 
ciation with a given 2N potential model. We report here our progress in the first task. Contact 
terms have to respect the discrete symmetries of the underlying theory, namely (besides parity) 
charge conjugation, or better to say, in a non-relativistic theory, time-reversal invariance. In addi- 
tion, strong constraints come from Fierz identities combined with the anticommutation properties 
of field operators, or equivalently, from the fact that the nuclear wave function is antisymmetric 
under exchange of two particles. As an illustration we consider in the next section the leading (no 
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derivatives) isospin-breaking contact operators. 

2. Leading isospin-breaking contact operators 



The most general isospin- violating, charge-invariant isospin structures of nucleons i,j,k are 

(2.1) 



T- = (x i xx j ) i , 

where the superscript refer to parity under time reversal. E ven (odd) isospin structures have to be 
associated with (rotation invariant) spin structures containing even (odd) numbers of a matrices. 
As a result, a list of seven operators is found 

0\ = N*T?NN*NN*N, (2.2) 

0\ = N^T 3 a 'NN^ a'NN^N , (2.3) 

0\ = NVAWVAWVAf, (2.4) 

0\ = e ab3 e ijk N t r a o i NN t T b o J NN t o k N, (2.5) 

0\ = nVawVawVn, (2.6) 

0\ = Nh"a i NN f z a a i NN i z 3 N, (2.7) 

O] = NU a o i NN 1 x a NNW 1 o i N. (2.8) 

Simultaneous exchanges of spin and isospin indices allow to write a set of linear relations among 
these operators. For example, from the exchange of nucleon fields 1-2, 

0\ = -\o\-\o\, (2.9) 

Ol = ~\o\ + \ol (2.10) 

0\ = -\o\-\o\-\ol (2.11) 
0\ = 20\-20\, (2.12) 



3^.1 



O] = ~^Ol + -O r 7 . (2.15) 

Two further set of relations are found by exchanging the other two pairs of nucleons. It turns 
out that six relations are linearly independent, which implies that only one operator survives, for 
example 0\. 

3. Subleading isospin-conserving contact operators 

The same procedure can be applied to all possible isospin-invariant 3./V contact operators with 
two spatial derivatives JT4|]. The initial list of 146 operators is reduced to 14 after imposing the 
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Fierz constraints, which in this case also involve the field derivatives. A further reduction can 
be operated invoking Poincare covariance, which can be implemented order by order in the low- 



energy expansion, as done in [|15|]. At the order we are working, this requirement amounts to select 
operators which do not involve the total momentum 1 . As a result, a list of 10 independent operators 
is found, and the 3N contact potential in momentum space can be written as 



y= I 

¥J¥* 



-£jkf - E 2 kj ti ■ tj - E 3 k}Oj ■ Gj - E A kj<Ji ■ Oyfj • Zj 

-E 5 (3k,- • ok; ■ Oj - kj) - E 6 (3k; • <7,k, • Gj - kf) )t, ■ tj 

+ ^E 7 ki x (Q,- - Qj) ■ (Gi + Gj) + l -E % ki x (Q, - Qj) • (ff,- + aj)Xj ■ X k 



-E 9 kj ■ Gikj ■ Gj - Eioki ■ Gikj ■ CTyT; • Ty 



(3.1) 



where k, = pi — pj and Q, = pi + pj are given in terms of the initial and final momenta of nucleon 
i, respectively p, and pj. Using a momentum cutoff depending only on momentum transfers, as 



done in [jl6p, e.g. F(kj; A)F(k|;A), the ensuing coordinate space 3./V potential has a local form, in 
terms of the function 



z »( r;A) =/w e ' prF(|,2;A) 



(3.2) 



and its derivatives. Explicitly, 



V= £ (E l +E2T i -Xj+E 3 G i -Gj + E 4 X i -i: j Gi-Gj 

+(E 5 +E 6 TrT j )S i j\z^r ij ) 



Z'o(rij)+2 



Zo(ra) 



Zo(r ik ) 



+(E 1 +E s T i -T k )(L-S) 



Z'(r, 



OViJ, 



Zo(r ik ) 



+ (E 9 +E w Zj ■ Z k )Gj ■ TijGk ■ r ik Z' (r u )Z' (r ik ) 



(3.3) 



where S;y and (L • S),y are respectively the tensor and spin-orbit operators for particles i and j. 
As it is apparent, a choice of basis has been made such that most terms in the potential can be 
viewed as an ordinaiy interaction of particles ij with a further dependence on the coordinate of 
the third particle. In particular, the terms proportional to Ej and £g are of spin-orbit type, and, as 



already suggested in the literature [17], have the right properties to solve the A y puzzle (see also 



Ref. [18, tt 



4. Matching to the pionful theory 

In the pionful theory the 3N interaction starts at N2LO, then it is parameter-free at N3LO 
[12, 13 ], and it's being currently investigated at N4LO [ 20 ] . In the very low-energy regime p M n 
pion exchange diagrams reduce to 3./V contact terms. This is achieved by Taylor-expanding the 



We ignore for the time being the "fixed terms", whose couplings are determined in terms of the lowest-order LECs. 
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corresponding 3N potential in momentum space. For example, the two-pion exchange component 
of the 3N potential at N2LO, reduce at low energy, to 

2 

V£ 2L0 ~ -2§§2 ff i kiff 2 k 2 Tj ■ T 2 + perm. + 0(k 4 ), (4.1) 
while the one-pion exchange component yields 

W^-con, ~ • k2ff2 • k2Tl • T 2 + perm. + 0(k 4 ). (4.2) 

Afterwards, the induced operators have to be expressed in the chosen basis, using the explicit rela- 



tions provided in [14]. An analogous procedure can be applied to the N3LO expressions provided 



in Refs. [12, 13], in a straightforward, although tedious manner. As a result, one gets the pion- 



induced contributions to the E\ \q as functions of the lowest order LECs, 

Ex = 755g ^ + 8 1 &± ^_~010 (43) 

2A516%F%M n 256nF%M n 64nF*M n 8mF%Af%. ' ' v ' ' 

E 2 = + 2 ^ Ct - 5 & Ct - ^ 5Ct + 2Cs K 006 (44) 

3686471^^ 384^ Mjr 19271^^ 48mF 2 M 2 ' ' v " 

E , = 3 ^a | 3glC 7 9g 2 C r 

2048wF|Af„ 64lcF*M* 16mF 2 M 2 

E ^ = -mS^-j^r ^ (4 - 6) 

79 ^ . + d 8 * Ct ^o,Q2, (4.7) 



12288^Af ff 256nF^M„ 6AnF%M K &mF%M% 
E 319g^ gj gl(Q-2C r ) 

6 36864^^ 256nF$M„ 24m7rF 2 M 2 ' ' V ' 

E 7 = ?M ^— + ~-0 08 (4 9) 

6\44kF%M k \2%nF%M n 4mF%M% ' ' 

^ 8 = ~ 3Xfl2nF$M n ~ \2%%F%M n + 4mF 2 M 2 ~ ~°-° 2 ' (410) 

19 ^ ^ - 0.14, (4.11) 



4096nF|M s 8mF 2 M 2 
c lg 2 g A D 427g^ 9gl gl(C 5 + C r ) 

10 2F*M 2 8F 2 M 2 12288^F„ 6 M ff 5\2nF$M x &mF%M% ■ ■ K ■ ) 

The (rough) numerical values are given in units of F%M„. Notice that only E\o receives contri- 
butions from both N2LO (~ —0.05) and N3LO (~ 0.15) and there is no sign of convergence. At 
N4LO there will appear the "genuine" contact contributions: if the convergence pattern is so bad, it 
is possible that they be phenomenologically relevant. In any case, the choice of a purely contact 3N 
interaction (including the subleading one), with a cutoff ~ M n , could turn out to be rather effective 
from the practical point of view. 
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5. Outlook 

In view of the above remarks, although in the chiral theory the subleading contact terms would 
only contribute at N4LO, they could play an important role, at sufficiently low energy, and provide 
the necessary flexibility to improve the description of the scattering observables, in particular the 
N — d vector analyzing powers A y and iT\\ and the tensor analyzing power T21, as well as several 



breakup cross sections Q21Q. Having fitted, in conjonction with a given realistic two-nucleon inter- 
action, these 3N observables, the single AN contact term [ p] ] could be used to reproduce the 4 He 
binding energy, and the so-determined interactions tested in four-body scattering states. 
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